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When a neutral atom moves in a properly designed laser field, its center-of-mass motion may 
mimic the dynamics of a charged particle in a magnetic field, with the emergence of a Lorentz-like 
force. In this Colloquium we present the physical principles at the basis of this artificial (synthetic) 
magnetism and relate the corresponding Aharonov-Bohm phase to the Berry's phase that emerges 
when the atom follows adiabatically one of the dressed states of the atom-laser interaction. We 
also discuss some manifestations of artificial magnetism for a cold quantum gas, in particular 
in terms of vortex nucleation. We then generalise our analysis to the simulation of non-Abelian 
gauge potentials and present some striking consequences, such as the emergence of an effective 
spin-orbit coupling. We address both the case of bulk gases and discrete systems, where atoms 
are trapped in an optical lattice. 
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Introduction 



In 1982 Feynman introduced the concept of a quan- 
tum emulator, as a possibility to circumvent the diffi- 
culty of simulating quantum physics with classical com- 
puters ( Feynman] |l 982 ) . His idea, based on the univer- 
sality of quantum mechanics, was to use one controllable 
device to simulate other systems of interest. Nowadays 
Feynman's intuition is being implemented in various se- 
tups and among them, cold gases of neutral atoms play a 
central role ( |Bloch et al \ |2008| |Buluta and Nori[ [2009 1 . 



* Electronic address: jean.dalibard@lkb.ens.fr 
^Electronic address: fabricc.gcrbicr@lkb.cns.fr 
tElectronic address: Gcdiminas.Juzcliunas@tfai.vu.lt 
§ Electronic address: P.Ohberg@ hw.ac.uk| 



These gases constitute remarkably flexible playgrounds. 
They can be formed of bosons, fermions, or mixtures of 
both. Their environment can be controlled using the po- 
tential created by laser light, with harmonic, periodic, 
quasi-periodic or disordered energy landscapes. Interac- 
tions between particles can be adjusted using scattering 
resonances. At first sight the only missing ingredient for 
implementing Feynman's idea with dilute atomic gases is 
the equivalent of orbital magnetism, which would allow 
one to simulate phenomena such as Quantum Hall effect. 
Here we discuss a general path that has recently emerged 
to fill this missing item, by using atom-light interaction to 
generate artificial gauge potentials acting on neutral mat- 



ter. The major part of this Colloquium addresses the case 
of a gauge field which has the U(l) Abelian symmetry, 
like electromagnetism, but we also discuss more elaborate 
configurations leading to non-Abelian gauge potentials. 

From the quantum mechanics point of view, the orbital 
magnetism of a particle with charge e can be viewed as 
a consequence of the Aharonov-Bohm phase 7 acquired 
by the particle when it travels along a closed contour C 
(Aharonov and Bohm 19591. This phase has a geomet- 



ric origin and does not depend on the duration needed 
for completing the trajectory. It reads 7 = 27r$/<I>o, 
where $ is the flux of the magnetic field through the 
contour C and $0 = h/e is the flux quantum. Therefore 
the quest for artificial magnetism amounts to realize sit- 
uations where a neutral particle acquires a geometrical 
phase when it follows the contour C. Similarly the gener- 
ation of a non-Abelian field can be achieved for particles 
with an internal degree of freedom. After completion of 
the trajectory along C, the initial internal state of the 
particle is changed into U\tj)i), where U is a unitary 
operator acting in the internal Hilbcrt space of the par- 
ticle and depending only on the geometry of the contour. 
Note that through all this Colloquium the artificial gauge 
potentials that appear are not dynamical variables, i.e. 
they are not influenced by the motion of the atoms. 

One of the simplest examples of geometric phases, the 
so-called Berry's phase, plays a central role in this pa- 
per. It emerges for example when a neutral particle with 
magnetic moment fj, moves in a (real) non-homogenous 
magnetic field Bq(t) (Berry 1984). If the particle is pre- 



pared at a point r in one of the local eigenstates \m(r )) 
of the Hamiltonian — fi ■ B (r) and moves slowly enough, 
it follows adiabatically the local eigenstate \m(rt)}. Once 
the trajectory along C is completed, the particle returns 
to the internal state |m(ro)), up to a phase factor con- 
taining a geometric component. 

Berry's phase is also present in atom-light interactions 
( |Dum and Qlshanii[ [1996). The role of the magnetic 
states \m(r)) is now played by the dressed states, i.e. 
the eigenstates of the atom-light coupling. The dressed 
states can vary on a short spatial scale (typically the 
wavelength of light) and the artificial gauge fields can be 
quite intense, in the sense that the geometric phase is 
large compared to 27r for a contour encircling a gas of 
realistic size. If the gas is superfluid, it will thus contain 
several vortices at equilibrium. The situation considered 
here should not be confused with phase-imprinting meth- 
ods, which correspond to imposing a dynamical rather 
than a geometrical phase to the atomic wavefunction [see 
Andersen et al. ( |2006[ ) and refs. in]. In the latter case 
one acts on the atoms with a time-dependent atom-light 
coupling, and the spatial phase profile of the light beam is 
dynamically transferred to the atom cloud. Here by con- 
trast we look for a time-independent Hamiltonian and 
require that its ground state has a non-zero vorticity. 

The use of geometrical phases is not the only way to 
reach a Hamiltonian with an artificial gauge field for 
neutral particles. Rotating the system at angular fre- 



quency Q, around a given axis (say z) also leads to ar- 
tificial magnetism in the rotating frame, with B z oc Q. 
This has been widely u sed i n the context of quantum 
gases [see |Cooper| ([2008) and |Fetter| ( [2009] ) for reviews 
and it is well suited when the confining potential is ro- 
tationally invariant around the z axis. The Hamilto- 
nian is then time-independent in the rotating frame, and 
the standard formalism of equilibrium statistical physics 
is applicable. However if the confining potential has a 
non-zero anisotropy in the lab frame, the Hamiltonian 
is time-dependent in any frame and the situation is dif- 
ficult to handle from a theoretical point of view. By 
contrast, methods based on geometrical phases do not 
impose any constraint on the symmetry properties of the 
initial Hamiltonian and the artificial gauge fields are pro- 
duced in the laboratory frame. Together with the possi- 
bility to extend the scheme to non-Abelian gauge poten- 
tials, this represents a significant advantage. Of course 
the use of laser beams also comes with some drawbacks 
that we will meet in this Colloquium, such as heating of 
the atoms because of residual spontaneous emission. 

This Colloquium is organized as follows. In Sec. [I] we 
present a toy model, where the atomic transition that 
is coupled to the laser light is represented by a simple 
two-level system. It allows us to identify the important 
elements for the emergence of artificial gauge potentials: 
gradients of the phase of the light beam and of its inten- 
sity or detuning with respect to atomic resonance. This 
model is relevant for some atomic species like Ytterbium 
atoms that possess an electronic excited state with an 
extremely long lifetime, but it cannot be used as such 
for the the alkali-metal atoms, which are the most com- 
monly used species in experiments. Indeed the large rate 
of spontaneous emission processes and the ensuing ran- 
dom recoils of the atoms cause in this case a prohibitive 
heating of the gas. In Sec. [H] we thus turn to schemes 
where the ground state manifold is degenerate, first in- 
troduced by Dum and Olshanii ( 1996), and we show that 



a gauge potential with a non-vanishing curl can then be 
generated even if the population of the e xcited states is 
negligible (Juzeliiinas and Ohberg 2004 1. We also dis- 



cuss the recent implementation of an artificial magnetic 
field in a Bose-Einstein condensate of rubidium atoms by 
Lin et al. (2009a|), which led to the observation of quan- 



tized vortices. Section III is devoted to the preparation 



of non-Abelian gauge fields and to the discussion of some 
of their physical consequences: negative refraction and 
reflection, implementation of the Klein paradox and the 



non-Abelian Aharonov-Bohm effect. In Sec. IV we study 
the production of artificial gauge potentials in optical lat- 
tices. Starting from the proposal by |Jaksch and Zoller] 
( 2003 ) , we show that one can reach the strong magnetic 



field regime, where the phase 7 per lattice cell can take 
any value between and 2ir. We also briefly discuss the 
implementation of non-Abelian schemes in a lattice. 

Throughout this Colloquium, we mostly address single 
atom physics, with some incursions to many-body physics 
at the mean-field level when dealing with the generation 
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I TOY MODEL 



of quantized vortices in a superfluid. The application of 
a gauge field on a quantum gas can also lead to the emer- 
gence of strongly correlated many-body states. Their de- 
tailed discussion is outside the scope of this Colloquium, 
and we simply indicate a few lines of research in this 
direction in the final Outlook section. 



I. TOY MODEL: TWO-LEVEL ATOM IN A LIGHT BEAM 



the Hamiltonian Q-Q with an alkaline-earth atom ir- 
radiated by a quasi-resonant laser beam, and we discuss 
the physical origin of the gauge potentials. We also study 
under which condition the strength and the spatial extent 
of the geometric magnetic field are sufficient to induce a 
large circulation of the atomic phase. This can then al- 
low for the nucleation of a lattice of quantized vortices if 
one applies this scheme to a collection of identical atoms 
forming a superfluid. 



In order to present the essential ingredients of the 
physics of geometrical gauge fields, we start this Col- 
loqium by discussing the simplest scheme for which arti- 
ficial magnetism can occur. We consider a single quan- 
tum particle with a two-level internal structure, and we 
show how the adiabatic following of one particular inter- 
nal state can provide the desired gauge fields. This will 
allow us to present with a minimum algebra the impor- 
tant physical concepts, which will subsequently be gen- 
eralized to more complex schemes. 

We denote {|g),|e)} a basis of the 2-dimensional 
Hilbert space associated with the internal degree of free- 
dom of the particle. Later on these states will represent 
electronically ground and excited states of an atom, re- 
spectively. We assume that the particle evolves in space- 
dependent external fields that couple \g) and |e). At the 
present stage we do not specify the physical origin of 
these fields. In practice these can be the fields of optical 
lasers, microwave fields and/or static electric or magnetic 
fields acting on the electric or magnetic dipole moment 
of the particle. The general form for the Hamiltonian of 
the particle of mass M is 



H = 



- (EE 

~ \2M 



V) 1 + U, 



(1) 



where P = — i?iV is the momentum operator and 1 is 
the identity operator in the internal Hilbert space. The 
coupling operator U can be written in the matrix form 



U 



mi 

~2~ 



cos o 
sin( 



(2) 



The particle dynamics is determined by four real quan- 
tities that may all depend on the position vector r: The 
potential V acts on the particle in a way that is indepen- 
dent of its internal state, whereas the generalized Rabi 
frequency f2 characterizes the strength of the coupling 
that lifts the degeneracy between \g) and |e). The two 
remaining quantities are the mixing angle 9 and the phase 
angle <fi. For a two-level atom in a monochromatic laser 



field [see Eq. (10) below], SI cos stands for the laser de- 
tuning from the atomic resonance, Jl sin 9 is the magni- 
tude of the atom-laser coupling and <p is the laser phase. 

In this section we first describe the atomic dynamics 
when the internal state of the particle follows adiabati- 
cally one of the eigenstates of U, and we give the expres- 
sion of the geometrical gauge potentials that appear in 
this case. Then we present a possible implementation of 



A. Adiabatic following of a dressed state 

At a point r the eigenstates of the matrix U are 



. , / cos(0/2) \ . 
|Xl}_ U*sin(0/2)J ' lX2 > 



i(0/2)\ 
cos(0/2) J ' 



(3) 

with eigenvalues Ml/2 and —Ml/2, respectively. We 
will call them dressed states, anticipating on the follow- 
ing discussions where they will correspond to the local 
eigenstates of the Hamiltonian describing the coupling 
between an atom and a light field. Since the states 
{ I Xi ) } form a normalized, orthogonal basis, the quantity 
i(Xj\^Xj) is a rea l number and the relation (VX2IX1) = 
-<X2|Vxi> holds, where we set \VXj) = V i\Xj))- 

Using the {|xj)} basis for the internal Hilbert space, 
we can write the full state-vector of the particle as 



|¥(r,i)>= ]T lfc(r,i)IXj(r)). 

j=l,2 



(4) 



Suppose now that at initial time the particle is prepared 
in one particular internal dressed state, say |xi). If the 
velocity distribution of the particle involves only suffi- 
ciently small components, the internal state of the parti- 
cle will remain proportional to \xi) for all times. This is 
equivalent to the Born-Oppenheimer approximation in 
molecular physics: The position r of the particle and its 
internal degree of freedom play the role of the nuclear 
coordinates and of the electron dynamics, respectively. 

We now derive the equation of motion for ipi in the 
case where "02 remains negligible at all time. We consider 
first the action of the momentum operator P on the full 
state-vector |*). Employing V [tpj \xj)] = [V^ 3 -] \xj) + 
tjjjl^Xj) an d the completeness relationship, one has 



P|*> 



(5) 



= ift(Xj\^Xl)- Assuming that ip 2 = 0, we 



with A 

project the Schrodinger equation iH\^f) = H\^f) onto the 
internal dressed state |xi)j where H is the full Hamilto- 
nian This leads to a closed equation for the proba- 
bility amplitude -01 to find the atom in the first dressed 



state (|Berry 


1989| Jackiw 1988; Mead 


1992 Mead and 


Truhlar[ 197S 


Moody et aZ.||1989|: 



ih 



chpi 
dt 



2M 



V 



~2~ 



w 



(6) 
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In addition to the terms V and Ml/ 2 that were already 
explicit in the initial Hamiltonian ( IJ , two geometric po- 
tentials A and W emerge in the adiabatic elimination of 
the state \X2), due to the position dependence of the in- 
ternal dressed states. The first one is the vector potential 

A(r)=tft(xi|Vxi> = |(coB0-l)V0. (7) 

The effective magnetic field associated with A is 

n 



B(r) = V x A = - V(cos6>) x Vc, 



(8) 



When this artificial magnetic field is nonzero, the vector 
potential A cannot be eliminated from Q by a gauge 
transformation. The particle acquires an effective charge 
(that we set equal to 1 by convention) and its motion 
exhibits the usual features associated with orbital mag- 
netism. A nonzero value of B occurs only if both the 
mixing angle 9 and the phase angle (f> vary in space with 
non collinear gradients. The second geometrical poten- 
tial appearing in ^ is the positive scalar term 



W(r) 



2M 



l(X2|Vxi 



8M 



[(V£ 



(9) 

The first experimental evidence for the scalar potential 
in Quantum Optics was given by Dutta et al. (1999). As 
such, the geometric scalar potential is not a very use- 
ful tool in the sense that there exist many other ways 
to create a scalar potential on atoms, using for example 
the AC Stark shift created by a far-detuned laser beam 



( Grimm et al.\ 2000 ) . However in the following we will es- 
timate the value of the scalar potential for some relevant 
configurations, because it must be taken into account to 
determine the equilibrium shape of the cloud. 

The contributions of the vector and scalar potentials 
m Eq. (jj) can be recovered in a systematic expansion 
in terms of the dimensionless adiabatic parameter, de- 
fined as the ratio between the short and long time scales 



(Littlejohn and Weigert 19931 Weigert and Littlejohn 



1993). The term —P ■ A/M appears at first order of the 



expansion, whereas A /2M and the scalar potential W 
appear at second order. An extra term that is quadratic 
with respect to P also shows up at second order, with a 
contribution that can reach that of the scalar potential 
when P 2 /M ~ Ml. We will not address this extra term 
in this review, because (i) we are mostly interested in 
the physics arising from the leading term of the expan- 
sion — P ■ A/M and (ii) we restrict to situations where 
the atomic kinetic energy is much smaller than Ml. 

The reason for which A, B and W are called 'geomet- 
rical' is clear from expressions 0, ([§ and which 
depend only on the spatial variation of the angles 6 and 
4>, i.e. of the geometry of the coupling between \g) and 
|e), but not on the strength fl of this coupling. Note 
that if we consider the adiabatic following of \\2) instead 
of |xi)j t ne equation of motion for ip 2 contains the same 
scalar potential W and the opposite vector potential A. 
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FIG. 1 (Color online) (a) A running wave propagating along 
the x axis, with a gaussian profile along the y axis (waist w) 
is used to create a geometrical gauge field on a two-level atom 
whose resonance frequency is close to the laser frequency, (b) 
Variation of the amplitude of the artificial magnetic field B 
(in units of Bq = hk/w) as a function of y/w; k'°'/A = 1 
(blue), 5 (red), 20 (black). 



B. Practical implementation with alkaline-earth atoms 

We now discuss how the above model can be imple- 
mented in quantum optics, in order to create artificial 



orbital magnetism on a gas of cold neutral atoms ( Dum 



and Olshanii 


1996 


Visser and Nienhuis 


1998) 



The sim- 
plest scheme consists in shining a single laser beam on 
an atom. We restrict the internal atomic dynamics to a 
two-level transition of frequency wa, between the ground 
state \g) and an electronically excited state |e). The laser 
light of frequency wl is supposed to be close to resonance 
with this transition. We suppose that the rate of spon- 
taneous emission of photons from the excited state |e) 
is negligible on the relevant time scale, which is a re- 
alistic assumption if the experiment is performed using 
the intercombination line of alkali-earth atoms (calcium, 
strontium) or Ytterbium atoms. Indeed for these atomic 
species, the radiative lifetime of the excited state |e) in- 
volved in the intercombination line [also used in optical 
atomic clocks ( Ye et al. 2008 )] is several seconds or even 



tens of seconds, large compared to the typical duration 
of cold atoms experiments. 

We suppose that the atomic center-of-mass motion is 
restricted to the xy plane using a suitable trapping po- 
tential that freezes the z degree of freedom. The mode of 
the laser beam is a Gaussian traveling wave of wavenum- 
ber k and wavelength A = 2n/k propagating along the 
x axis, with a waist w in the y direction (Fig. [T^.). The 
states \g) and \e) stand for eigenstates of the internal 
atomic Hamiltonian in the absence of coupling with the 
radiation field. Using the rotating wave approximation, 
the coupling matrix U can be written in the present case 



(Cohen-Tannoudji et al. 1992) 



U 



h ( A k* 

2 \ k -A 



(10) 



where the detuning A = wl — and where the Rabi 
frequency k characterizes the strength of the atom-laser 
coupling. The spatially varying phase of the laser is in- 
corporated into k which is therefore complex. 
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In the following we neglect the diffraction of the laser 
beam and set «(r) = k(y) e lkx , where k is real and pos- 
itive. The angle (j> featured in Eq. ^ is then simply 
the running phase kx of the propagating laser beam and 
V</> = ke x , where e x is the unit vector in the x direc- 
tion. The mixing angle 9 is given by tan# = k/A and 
two options are available to provide a non-zero V#. One 
can either use a gradient of the Rabi frequency k pro- 
duced by a spatial variation in the laser intensity, or a 
gradient of the detuning A. In this subsection we discuss 
these two configurations, which we refer to as "grad-K" 
and "grad-A" configurations, respectively. 



a. Grad-K configuration. A possible scheme leading to 
this configuration is shown in Fig. [l] where we take ad- 
vantage of the transverse Gaussian profile of the laser 
beam. The Rabi frequency is k(y) = e~ v ' w , which 
leads to an effective magnetic field B that is parallel to 
the z axis, with the amplitude obtained using Eq. (18]): 



r-2 



B = B n 



Ak 

1W 



(11) 



with SI 2 = A 2 + k 2 . Here we set Bq = Hk/w, which gives 
the typical scale for the effective magnetic field appearing 
in this case. The characteristic length scale associated to 
orbital magnetism is the so-called magnetic length £g = 
(h/B) 1 / 2 , which gives the size of the elementary quantum 
cyclotron orbit. We obtain in this context £ B ~ (wA) 1 / 2 
(hence £ B > A) for B ~ B - 

The variation of B with y/w is plotted in Fig. [TJd for 
various values of the ratio n^/A. When k'°) / A ^s> 1 
the maximal value of B is obtained approximately at the 
point where k = A, i.e. 2/ max ~ w [log(K ( ^ - ) / A)] 1 / 2 w 
and J5 max ~ B y max j '(2y2w). At first sight the limit- 
ing regime k^/A — > oo is interesting since it leads to 
B mM — > oo. However the interval Ay over which B takes 
significant values tends to zero as 1/-B max , and eventually 

1/2 

becomes smaller than the magnetic length £ B oc 1/B max , 
making this regime 3> A not useful in practice. More 
generally we find that 



+oo 



Bdy 



hk 
Y 



v/l + ( K (o)/A) 2 



< 



hk 



(12) 



Consider for example the case k( )/A = 5. The maxi- 
mum value -Bmax ~ 0.45 B is obtained for \y\/w w 1.2, 
and B > Bo /4 over an interval of width ~ w/2. 
The calculation of the scalar potential leads to 



W(y) 



Er, 



« 2 (y) 



4 A 2 + k 2 (y) 



Ay 2 



A^ 



w A k 2 A 2 + k 2 (y) 



, (13) 



where we define the recoil energy Er = h 2 k 2 /2M as the 
kinetic energy of an atom initially at rest when it absorbs 
or emits a single photon. Realistic values of the waist 
w are such that kw ^> 1. When < a few A, the 



main contribution to the scalar potential (J9J) is the (V0) 2 
term, which corresponds to the first term in the bracket 
of Eq. ( 13 ). The scalar potential then creates a potential 



bump that is maximal in y = 0, with a height scaling like 
_Er. When / A ^> 1, the contribution proportional to 
(V#) 2 to the scalar potential becomes dominant [second 
term in the bracket of Eq. ( |13| ]. The corresponding force 
is now very large around the point y max and may create 
strong distortions of the trapped atom cloud. 



b. Grad-A conRguration. We suppose now that the laser 
waist w is very large, so that the spatial mode of the 
laser is well approximated by a plane wave and the Rabi 
frequency k can be taken uniform. We also assume that 
a gradient of detuning is present along the y axis, so that 
A = A'(y — y ). For alkaline-earth or Ytterbium atoms, 
it can be created using an additional off-resonant laser, 
which produces a differential light shift for the g and e 
states. The effective magnetic field B is parallel to the z 
axis, with an amplitude given by 



B = B Q C 3 / 2 (y) 1 C(y) 



1 



i + (y - yo) 2 /il 



(14) 



where we introduced the characteristic length £ K — 
k/\A'\ and set Bq = hk/(2£ K ). The scalar potential is 



W(y) 



~4~ 



1 



fc 2 ^ 2 



C 2 (y) 



(15) 



The magnetic field can be made arbitrarily large by tak- 
ing A' — > oo , but as for the grad-K case this is not relevant 
in practice. Indeed the field then takes a large value only 
over a domain whose size £ K oc 1/A' becomes smaller 
than the magnetic length £ B = (^A/2) 1 / 2 . 



C. Validity of the adiabatic approximation 

We now discuss the validity of the adiabatic approxi- 
mation ( Messiah 1961b ) underlying the projection mech- 

For simplicity, 



anism generating the gauge potentials, 
we adopt a semi-classical point of view in which we treat 
classically the atomic center-of-mass motion, but still use 
quantum mechanics to describe its internal dynamics. 
The atom is supposed to be initially at rest in the inter- 
nal state |xi), and then accelerated in order to reach the 
velocity v. The population IJ 2 of the state |%2) at the end 
of the acceleration process is not strictly zero. For small 
\v\ it is given by IJ 2 ~ \v ■ (% 2 | Vxi)/^| 2 , where the mo- 
tional coupling is (xg lVxi) = [V0 - isin6» V^] e _i */2 
2008] ). 



(Cheneau et al. 



We first address the grad-K implementation discussed 



above, for which VO = AVk/(A 2 



J ), and we re- 



strict to the practical case where < a few A. Since 
the gradient of the laser intensity always occurs with 
the spatial scale w ^> we find that |V0| -C k. 

The motional coupling is dominated by the contribu- 
tion of and reaches the maximum value |(x2|Vxi) ~ 
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k\ sin(0)| = k\n(°) Using the fact that the atomic ve- 
locity changes by an amount equal to the recoil velocity 
vr = hk/M during each photon absorption or emission 
process, a necessary condition for the adiabatic approxi- 
mation to hold (1T2 <C 1) is thus 



In most situations of practical interest, 



Eq. ( 16 ) reduces to the intuitive condition 



/ik (0) > E R , 



(16) 
A and 

(17) 



stating that the coupling strength should be much larger 
than the recoil energy. Of course the validity condition 
for the adiabatic approximation is more stringent if the 
atomic velocity is notably larger than «r. 

We now address the case of the grad-A scheme. Far 
away from any resonance point where A vanishes, the 
adiabaticity condition ( 16 1 still holds. Near a resonance 
point, the validity condition for the adiabatic approxi- 
mation may be more stringent, since the gradient of the 
mixing angle |V0| ~ l/t K can be large. More precisely 
suppose that on each side of the resonance point, |A| is 
large compared to the coupling k. The dressed states 
then nearly coincide with the bare states, with cither 
{|Xi> = \e), IX2) = Iff)} or {| X i> = \g), |X2> = |e>}. The 
switch between these two configurations occurs within a 
region of size ~ l K around the resonance point. Taking 
again the atomic velocity on the order of the recoil ve- 
locity vr, the adiabatic condition reads 



hk- 



1 + (ki K y 



1/2 



(18) 



When kl K » 1, this reduces to Eq. (16) in y = y , i.e. 



hk 3> -Er. In the opposite limit kl K <C 1, the peak value 
of B is very large, but it is reached only in a small re- 
gion of width £ K <C fc" 1 and the adiabaticity condition 
kl K 3> wr is more demanding. This is thus not a conve- 
nient configuration for a continuum implementation, but 
it becomes relevant for the lattice case (Sec. IV). 



D. Physical interpretation of the geometric potentials 

We now turn to the physical interpretation of the ge- 
ometrical gauge fields. We focus our discussion on the 
two-level system of Sec. |I.A| and Sec. |I.B| but the follow- 
ing physical images can be generalized to the schemes 
that will be analysed in the next sections. 

The scalar potential W can be interpreted as the ki- 
netic energy associated with the fast micromotion of the 
particle. This was first explained for a classical contin- 
uous internal degree of freedom by |Aharonov and Stern 
( 1992 ) . Here we outline the reasoning of Cheneau et al. 



(2008) that addresses the case of a quantized internal de- 
gree of freedom. Consider a particle prepared in the in- 
ternal state \xi)i with a center-of-mass state that consists 



in a wave packet localized around a given point r with 
an extension small compared to the scale of variation of 
8 and 4>. We introduce the force operator F — — VC/ 
and note that the eigenstates of the coupling U are 
not eigenstates of F. The force acting on the particle 
thus exhibits quantum fluctuations, (-F 2 ) 7^ (F) 2 , which 
are characterized in the Heisenberg picture by the sym- 
metrized correlation function 



C(t) = ~(6F(0)-6F(t) + 6F(t)-6F(0)), 
= h 2 n 2 \( X 2\V X i)\ 2 cos(ftr). 



(19) 



with SF(t) = F(t) — (F) . We now recall that in classical 
mechanics, a particle submitted to a rapidly oscillating 
force 5F undergoes a micromotion with the average ki- 
netic energy 



E K = I ^du, 



(20) 



where C(w) is the Fourier transform of C(r). Inserting 
the value ( 19 ) of C(r) into ( 20 ) , we find that E^ coincides 
with the scalar geometric potential given in (p}]). 

The vector potential A given in Q is related to the 
Berry's phase 7 that appears when a quantum system 
- here, the two-state system associated with the inter- 
nal degree of freedom of the particle - is slowly trans- 
ported round a contour C, while remaining in an eigen- 
state \x( r )} of its Hamiltonian (Berry 1984 1 : 



7(C) = i j> (x| Vx> -dr=^<j)A-dr, 



(21) 



where the second equality holds for the two-level system 
considered above when it is prepared in the state \x}- 

When B = V x A ^ 0, the Lorentz force F = v x B 
acts on the atom when it moves with velocity v. For 



the scheme discussed in Sec. LB the momentum change 
imparted by the Lorentz force has a simple physical in- 
terpretation. Suppose that the atom moves along the 
y-axis, with a trajectory starting from y\ 3> w at time 
t\, and ending in yi — at time ti. Since B is directed 
along the z direction and v along the y direction, the 
average momentum change (AP) is directed along x: 



(AP, 



v y (t)Bdt 



ti 



Bdy. 



(22) 



Using (12) and y\ w, we find that AP X « hk/2 in the 
case where ^> A. This result has a simple physical 
interpretation. When the atom is located in yi 3> w, 
the occupied dressed state is \g). When the atom ar- 
rives at 2/2 = where the atom-laser coupling is maxi- 
mal, the dressed state is w (\g) + \e)e lkx ) /V2 in the limit 
3> A. Hence a measurement of the atomic momen- 
tum AP X can give the results or Hk, both with probabil- 
ity 1/2, hence (AP X ) — hk/2. A more detailed discussion 
of the physical origin of the Lorentz force for other atomic 
trajectories is given by Cheneau et al. (2008). 
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E. Achieving states with a nonzero circulation 

After this discussion of some simple schemes that gen- 
erate artificial gauge potentials, we now turn to the effect 
of this potential on the external (center-of-mass) degree 
of freedom of the atom. One of the main motivations 
for the generation of artificial magnetic fields is indeed 
to create some extended regions where the orbital mag- 
netism is sufficient to favor states with a nonzero circu- 
lation. For instance, if a superfluid is placed in such a 
region, its ground state will exhibit a vortex lattice. We 
now explore under which condition this can occur for the 
simple schemes outlined above. 

When a charged superfluid (with e = 1 here by con- 
vention) is placed in a magnetic field, the vortex density 
is p v = B/(2nh) or i n other words p v = i^ 2 , wher £g is 
the magnetic length ( [Tinkham 1996). If the magnetic 
field obtained from the geometric potential A keeps a 
value ~ B on a disk of radius r, one therefore expects 



that N v 



nr 2 p v 



B/(2H) vortices will be present 



in steady state in a superfluid filling this disk. We now 
wish to determine whether one can reach a situation with 
N v 3> 1, which is equivalent to requiring that the phase 



7(C) defined in (21 1 is large compared to 2n. 

Consider again the grad-K scheme represented in Fig.[T] 
where we choose for example K^- ' = 5 A. One gets in 
this way a fictitious magnetic field that is approximately 
uniform with a value ~ B Q /4 in the band parallel to the 
x axis, centered on y — 1.2 w, with a width £ y w w/2. 
The length of this band along x is limited only by the 
diffraction of the laser beam, which occurs on a distance 
3> w, if the waist w is chosen much larger than the laser 
wavelength A = 2ir/k. In order to study the physics of 
vortex lattices in this geometry, the requirement is thus 
simply to fit several vortex rows in the width £ y . Since 

— 1/2 

the distance between two vortex rows is ~ p v , this 
requirement can be written 



vortex rows 



Clearly this method is well adapted to the study of vor- 
tex lattices only if w ^> A. The choice of a small waist 
(of the order of a few A only) is not appropriate, because 
the scaling Bq oc 1/w, which would tend to favor small 
waists, is compensated by the other scaling £ y ~ w over 
which the field keeps a significant value. A similar argu- 
ment can be made for the grad-A scheme with £ y — £ K . 

So far we have restricted our discussion to the case of 
a single laser travelling wave and the spatial scale of vari- 
ation for the mixing angle 8 is thus the beam waist w. 
It is interesting to consider also the case where several 
traveling waves irradiate the atom at different angles, 
so that interference phenomena can introduce a much 
shorter length scale for 9, typically A/(27r). For simplic- 
ity we restrict to the case of two waves and we choose the 
corresponding wavevectors equal to k± — k(e x ±e y )/\/2. 
The resulting light field still has a spatially varying 
phase <f) — kx/\/2, and it presents an interference pat- 



tern along the y direction with a spatial period X/y/2. 
Hence |V>| ~ |V0| ~ k/\/2, and we find using Q that 
the maximal modulus of the artificial magnetic field is 
\B\ ~ 0.1 hk 2 K^/ A. This field is directed along the 
z axis, and is a periodic function of y with changes of 
sign every \/{2^/2). The same reasoning as above shows 
that one can marginally localize one quantum of circu- 
lation in each disk of area fc -2 over which the field is 
approximately uniform. In order to obtain a circulation 
2ir, one needs to rectify this spatially alternating field. 
Practical solutions will be detailed in Sec. IIVI devoted to 
artificial gauge fields in optical lattices. 



II. GAUGE POTENTIALS FOR MULTI-LEVEL SYSTEMS 

In the model discussed in Sec. [ij the internal state of 
the atom is at any place a linear combination of ground 
and excited states, and each of these two states must 
have a relatively large weight in order to obtain a non 
negligible artificial gauge potential. Therefore this con- 
figuration can be used only if the excited electronic state 
has a very long radiative lifetime, as is the case for 
alkaline-earth species. In order to address a larger class of 
atoms (including the more widely used alkali atoms), we 
now turn to schemes that take advantage of the (quasi- 
) degeneracy of the electronic ground level. Denoting 
{\dj)i j = 1, - ■ ■ , ^V} a basis set of the ground state 
manifold, we look for configurations where some dressed 
states are linear combinations of the \gj) states, with 
a negligible contribution of the excited state manifolds, 
|x) ~ Y^j a j\9j)- A s we wm see ' tms can be obtained by 



taking benefit of a so-called dark state ( Arimondol 1996 1, 



or by choosing a laser frequency that is strongly detuned 
with respect to the atomic resonance lines. If the atom 
is prepared in such a dressed state and moves sufficiently 
slowly to follow it adiaba tical l y, geometrical gauge poten- 
tials show up as in Sec. I.B (Dum and Olshanii, 1996 1 . 



(23) Since we use laser beams to provide the relevant stim 
ulated Raman couplings between the states \gj), the ctj 



coefficients can vary significantly on a short length scale, 
typically an optical wavelength. One can thus produce 
geometrical fields with comparable amplitudes to those 
found in Sec. |TJ while avoiding the strong heating that 
would be caused by spontaneous emission processes. 

In this section we first consider the dark state case, 
which occurs for a A level scheme, where two sublevels of 
the electronic ground state \gi) and I52) are coupled to a 
single excited state |e) by two laser beams. The dark 
state is an eigenstate of the atom-laser coupling that 
is a linear combination of \gi) and I32) with a strictly 
zero contribution of the excited state. We then discuss 
two possible practical implementations of this dark state 
scheme, first using laser beams carrying orbital angular 
momentum, and then using counterpropagating Gaus- 
sian beams with a spatial shift of their axis. Finally 
we describe an alternative scheme involving a position- 
dependent detuning. This scheme that is not relying on 



7 



U MULTI-LEVEL SYSTEMS 




FIG. 2 Atomic A-level structure providing a dark state that 
depends parametrically on the Rabi frequencies «i and K2- 



dark states but on a large detuning, has led to the first 
experimental observation by Lin et al. (2009a) of a geo- 



metric magnetic field in the context of cold atom physics. 



A. Artificial magnetic field in a A scheme 

We consider the A-type atomic level structure repre- 
sented in Fig. [2j where two laser beams couple the atomic 
states | <7i ) and \g 2 ) to the third one |e). The lasers are 
tuned symmetric with respect to the average of the fre- 
quencies of the gi — e and 52 — e transitions. The full 
atomic Hamiltonian is given in Eq. ([I]), and the coupling 
operator between the light and the atom written in the 
{Iffi): \ e )^ 1 92)} basis reads using the rotating wave ap- 
proximation 



U = 




(24) 



Here K\ t 2 are the complex, space-dependent Rabi fre- 
quencies, which include the spatially varying phases of 
the laser beams as in Sec. |I.B| The frequency 28 is the 
detuning of the two-photon excitation with respect to the 
Raman resonance between g\ and g 2 - 

Suppose that the two-photon (Raman) excitation is 
resonant (6 = 0). In this case the coupling matrix U 
possesses an eigenstate with zero energy called dark (or 
uncoupled). This state contains no contribution from the 
excited state |e) and reads: 



\D) = (k 2 |.9i) - ki|32))/«, 



(25) 



where k — (|ki| 2 + \k 2 



The two other eigenstates 



of U have the eigenenergies ±hn/2, and read |±) = (\B)± 
|e))/\/2, where \B) is the bright (coupled) state 



\B) = K\gi) + K* 2 \g 2 ))/K. 



(26) 



Dark states are frequently encountered in quantum op- 



tics a pplications such as sub-recoil cooling ( |Aspect et al. 
1988), electromagnetically induced transparency (Ari- 



mondo 



Lukin 



1996 Flcischhaucr et al. 2005 Harris 1997 



2003 ) and Stimulated Raman Adiabatic Passage 



state \D) with respect to the decoherence caused by spon- 
taneous emission. 

Like in Sec. |TJ the full atomic state-vector can be cast 
into the eigenstates of the operator U as 



|*(r)) 



X=D,± 



ip x (r)\X(r)} , 



(27) 



where the wave- functions V , d( t ") an d r( P±{ r ) describe the 
translational motion of an atom in the internal states 
\D(r)) and | ± (r)), respectively. The adiabatic approxi- 
mation assumes that the atom stays in the dark state, so 
that one can write approximately |^(r)) « ipr)( r )\D(r)) . 
Projecting the Schrodinger equation onto the dark state 
and neglecting the couplings to the two other internal 
states |±), we arrive at an equation of motion for the 
dark state wave function ipD^) similar to ([6]): 



(P-A) 
2M 



V + W 



(28) 



where A = ih(D\VD) and W = h 2 \(B\VD)\ 2 j(2M) are 
the effective vector and scalar potentials emerging due to 
the spatial dependence of the dark state. 

We therefore recover a situation similar to that of 
Sec. LA provided we replace \x\) by \D) and \x%) by 



B). Comparing Eq. ^ and Eq. (25), we set: 



7C: 



- — r = — tan - 

M 2 



(29) 



where the (f>j are the phases of the Rabi frequencies Kj — 
,j (j = 1,2). The artificial magnetic field B = V x 
A given in pj) can be expressed in terms of ( and <p: 



oi<t>. 



B 



_V0 x vc 
l (i + C) 2 



(30) 



This effective magnetic field B is non-zero only if the 
gradients of the intensity ratio £ and the relative phase 
4> are both non-zero and not parallel to each other. We 
discuss in the next subsections some practical implemen- 
tations of this A scheme, using either light beams with 
orbital angular momentum or counterpropagating Gaus- 
sian beams with an axis offset. 



B. Using light beams with orbital angular momentum 

We consider here the situation where the atoms lo- 
cated in the plane z = are irradiated by two laser 
beams propagating along the z axis (Fig. [3]) . The beams 
are prepared in Laguerre-Gauss modes, and they carry 
the orbital angular momenta and Mr> per photon. 
This scheme was first proposed by 



(|Bergmann et al 



2001). These app 



1998 Krai et al. 2007 Vitanov et al. 



Juzeliunas and Ohberg 
(2004) [see also (Juzeliunas et al. 2005a|b Zhang et al. 
2005)]. The complex Rabi frequencies can be written 



r) = Kj(p) e lljV (j = 1,2), where ip is the azimuthal 
ications rely on the robustness of the angle around the z axis and p is the radial coordinate 
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«2 




z 



Gas of ultracold atoms 



FIG. 3 (Color online) Two copropagating Gaussian beams, 
one of them being prepared in a Laguerre-Gauss mode with a 
non-zero orbital angular momentum, drive the two transitions 
of an atom with a A level scheme. For an atom prepared in 
the dark state \D) (Eq. 25 1, the nontrivial phase and intensity 
ratios between the Rabi frequencies Ki and k,2 produce an 
artificial magnetic field parallel to the beam propagation axis. 



in the x — y plane. The effective magnetic field (|30|) is 
z direct 

B(p) = 



directed along the z direction and its amplitude reads 

M d p C 



p (i + C) 2 



(31) 



where £ = l\ — £2 is the relative winding number of the 
two beams. 

For concreteness we consider beams with equal waists 
w and we choose £\ — £, £2 — so that 

K 1 (p) = KMp £ e-» 2 ' w2 , k2{p) = n^pie-P 2 / w \ (32) 

where p c is the radius at which the two beams have equal 
intensities. The winding number £ determines the shape 
of the gauge potentials, (i) For £ > 1, the effective mag- 
netic field is zero at the center of the beams and is max- 
imum at p = p c [{£ - I) J {I + l)] 1/21 . (ii) For t = 1, the 
effective magnetic field is maximum at the origin with 
the value B(0) = 2h/ p 2 . Its amplitude is inversely pro- 
portional to p c , and can thus be controlled by chang- 
ing the intensity ratio of the laser beams. Typically p c 
has the same order of magnitude as the waist w. Since 
w ^> fc _1 , the magnetic field generated in this way is no- 
tably smaller than the value Bq ~ hk/w found in Sec. LB 



where we took advantage of the rapid spatial variation of 
the phase kx of a plane running wave. 

The effective magnetic flux (in units of the Planck con- 
stant h) through a circle C of radius tq reads [see Eq. ( 21 1] 



7(C) 
2tt 



A ■ dl = £ ■ 



Co 



1 + Co 



(33) 



where Co is the intensity ratio at the radius p — ro- As 
explained in Sec. |I.E[ this flux gives the maximal number 
of vortices that can be observed if a superfluid gas is 
placed in this laser configuration. Its maximal value is £ 
and it is approximately reached for £0 3> 1, i.e. for radii 
Po such that the intensity of the beam 1 (with orbital 
angular momentum) largely exceeds that of the beam 2 
(with no angular momentum). In practice the winding 




FIG. 4 (Color online) Two counterpropagating beams with 
an offset 2a between their propagation axes drive the two 
transitions of an atom with a A level scheme. When the atom 
prepared in the dark state \D), the laser beams induce an 
artificial magnetic field perpendicular to the figure plane. 



number £ can reach a few tens (Hadzibabic 2011); this 



configuration is therefore better suited to generate small 
vortex patterns rather than large vortex arrays. 



C. Using spatially shifted laser beams 

We now turn to the scheme represented in Fig.[4j which 
constitutes a direct generalisation of the configuration 
studied in section II. Bl for a two-level atom. It uses two 
Gaussian beams that are counterpropagating along the x 
axis, so that the phase difference between the two beams 
provides the necessary gradient of the phase angle cj> en- 
tering in Eq. (30). The gradient of the intensity ratio C 



is obtained by a spatial shift ±a of each beam axis along 
the y direction. This configuration was first proposed 
by Juzeliunas et al. (2006), and a modified version that 



is more flexible in terms of choice of laser polarizations 
was later suggested by Giinter et al. (2009). It offers 



the possibility to reach large magnetic field values over 
an extended region of space, while taking advantage of 
the dark state configuration to minimize the spontaneous 
emission rate. This configuration with spatially shifted 
laser beams can also be used to generate the spin-Hall 
effect (Zhu et all 2006) and the Stern-Gerlach effect for 



chiral molecules (Li et al. 2007) 



To simplify the mathematical treatment we assume 
that both beams have equal waists w and equal cen- 
tral Rabi frequency k^°\ and we neglect their diffraction 
along the x axis. Furthermore we use the fact that the 
frequencies of the transitions g\ — e and g2 — e are very 
close in practice and we set k\ — &2 = k. The complex 
Rabi frequencies are then given by 



Kj {r) 



_ R (0) e ±ikx e -(y± a y/w 2 



(34) 



where the + (resp. — ) sign stands for j — 1 (resp j — 
2). The relative phase between the two beams is 4> = 
2kx, and the intensity ratio reads C = exp(8ya/w 2 ). The 
effective magnetic field is oriented along the z direction 
and its amplitude is obtained from Eq. |30|): 



B(y) 



Ahkc 



1 



w 2 cosh 2 (4ya/w 2 ) 



(35) 



The offset a of the beams along the y axis is a priori ar- 
bitrary. In practice, one should not take s > w, in order 
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excited state manifold 




nip = — 1 nip = 



trip = +1 



Lin et al. 



(2009a). 8 'Rb atoms 



FIG. 5 Atomic scheme used in 
with a spin F = 1 ground state are placed in an external mag- 
netic field that displaces the sublevels \m = ±1} (dotted lines) 
with respect to \m — 0). The atoms are irradiated by two 
laser beams of wave vectors fci and &2 that induce a resonant 
Raman coupling Amj = ±1 between the Zeeman sublevels. 
An additional spatial gradient of the external magnetic field 
induces an additional displacement of the \m = ±1} sublevels 
(full line). It thus creates a spatially varying two-photon de- 
tuning S which allows one to generate an artificial gauge po- 
tential. Note that the frequency scale is not respected in the 
drawing (in practice A 5). 



to keep a significant laser intensity along the line y = 
where B is maximum, and ensure that the atoms follow 
adiabatically the dark state \D) at this location. Taking 
as a typical value a — w/2 we find B(0) — 2Kk/w, which 
is comparable to the result of Sec. |I.B| Therefore the con- 
clusions of Sec. |I.E| concerning the possibility to generate 
large vortex arrays remain valid for this configuration. 
Like the geometric magnetic field B(y), the scalar poten- 
tial W(y) is maximum along the line y = 0. An extra 
trapping potential V is therefore needed to prevent the 
atoms from flying away from this area. 



D. Gauge potentials involving a gradient of detuning 

We explained in Sec. |I.B| that the necessary gradient of 
the mixing angle 9 can be provided by a gradient of the 
detuning of the laser frequency, as well as by a gradient 
of the laser intensity. The same distinction applies to the 
A scheme and more generally to schemes involving mul- 
tiple atomic levels. The first observation of a geometrical 
magnetic field by Lin et al. ( |2009a[ ) was actually based 
on a gradient of detuning for optical Raman transitions 
occurring between various ground sublevels. In this sec- 
tion we present the main features of this experiment and 
connect it with the already discussed configurations. 

The experiment by Lin et al. (2009a) has been per- 
formed with 87 Rb atoms in their F = 1 hyperfine level. 
In the process all three Zeeman sublevels \mp) with 
mp = 0,±1 acquire a significant population. The atoms 



are irradiated by two laser beams with wave vectors k\ 
and &2 that create a quasi-resonant Raman coupling be- 
tween Zeeman sublevels with Amp = ±1 (see Fig. [5j 
and its caption for details). The coupling occurs via the 
absorption of a photon in one beam and the stimulated 
emission of a photon in the other beam, accompanied by 
a change of the atom momentum of ±hk^, where the dif- 
ference of wave vectors fed = fci — ki = k ( ie x is chosen 
along x. An important ingredient is the application of a 
real magnetic field in addition to the laser beams. The 
Zeeman shift created by the magnetic field lifts the ini- 
tial degeneracy of the sublevels mp — 0, ±1. It gives a 
control knob on the two-photon detuning S. In the initial 
experiment by Lin et al. ( 2009b ) a spatially uniform real 



magnetic field was applied, and it led to a spatially uni- 
form vector potential, corresponding to a zero geometric 
magnetic field. Subsequently Lin et al. (2009a) used a 



non-homogeneous real magnetic field, making the two- 
photon detuning 6 position-dependent. In the following 
we assume the linear form 5 = S'(y — yo), where §' > is 
the uniform gradient of detuning. 

This experiment can first be analyzed using the adia- 
batic framework used so far. We assume that the cou- 
pling lasers can be considered as plane waves on the scale 
of the atomic cloud, so that the problem is translationally 
invariant along x. Therefore the two-photon coupling k 
induced by the laser pair can be written as k 



with <j) — fed • t = k^x being the phase of the Raman cou- 
pling. The single-photon detuning A with respect to the 
excited state manifold involved in the transition is sup- 
posed to be very large compared to the Rabi frequencies. 
One can then perform an adiabatic elimination of the ex- 
cited states to obtain an effective Hamiltonian acting on 
the ground state manifold. Keeping only the terms rele- 
vant for the subsequent discussion, this coupling written 
in the basis {\mp = —1), \mp = 0), \mp — +1)} has the 



form of Eq. ( 24 ) with ft 2 



We focus on the eigenstate of U associated to the lowest 
eigenvalue -h[6 2 + (ft(°)) 2 /2] 1/2 , 



|x)=e**coB a s|-l) 



^lO+e-^sin^l + l) (36) 



where we have set tan6> = k,(°'/v25. The vector poten- 
tial is now given by A = ih{x\^X) — —fifed cos The 
artificial magnetic field can be written in a form that is 



reminiscent of the result for the grad-A case of Sec. LB 
[Eq. (14)] (subject to the replacement of A by 8): 



B = e z B £ 3 / 2 (y-y Q ), 



(37) 



where the characteristic length entering in the Lorentzian 
function C is l K — /t' ' /V2S' ', and where the peak value 
of the artificial magnetic field is B Q = hk&jl K . The pre- 
diction (371 is in good agreement with the experimental 
results of |Lin et al. (2009a) {i K s» 40 A in the experi- 



ment). The scalar potential also takes a form that is 
similar to the one found for a two-level atom in the grad- 
A case: W(y) ~ fi. 2 fe^/(4M) C(y—yo), where we assumed 
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kd^K 3> 1- For an external trapping potential centered 
in yo, W(y) weakens the confinement along the y direc- 



tion by creating in the vicinity of yo the anti-trapping 
potential —Mu)y t/ {y — yo) 2 /%, with l o w = hkd/V2M£ K 
{low /27r ~ 30 Hz for the parameters of Lin et a/.| ( |2009a[ )). 

Alternatively, the scheme described above can be an- 
alyzed using an original framework introduced by | Spiel- 



man 



(2009). This framework is usable in the particu- 
lar case where the only x-dependence of the atom-laser 
coupling is the phase term e ±tkdX , and it has then a 
larger range of validity than the standard adiabatic ap- 
proximation. In a first step one "freezes" the y motion 
and diagonalizes exactly the x-dependent Hamiltonian 
H x = P*/(2M) + U , where U is the 3 x 3 matrix 
given in Eq. (24) with k 2 = «* = n^e lkdX . The eigen- 



states are 3-component spinors that can be labelled by 
their momentum KK X . For each K x , this diagonalization 
yields three dispersion relations E n (K x ) (n — 1,2,3), 
each with a minimum at a momentum K nym \ n depending 
on the two-photon detuning S and Rabi frequency n(°\ 
At low energy and large coupling strength one obtains 
E n s» h 2 (K x — K n>m i u ) 2 /2M* , where M* is an effective 
mass. In a second step, one considers the motion along 
y and takes into account the spatial variation of the de- 
tuning S — S'(y — yo). For the lowest energy state n = 1, 
this motion is described in an approximate manner by 
the projected Hamiltonian H = P*/(2M) + E\, which 
leads to the identification of an effective vector poten- 
tial A = hKi^ m i n e x . The spatial dependence of 6, hence 
of -KTi.min, along the y axis corresponds to an effective 
magnetic field B = H(dKi^ m i D /dy) e z , described here in 
the Landau gauge. The exact diagonalization and the 
adiabatic approaches agree when the latter is valid (see 
Sec. I.Cl, with M* « M in particular. For Rabi fre- 



quencies smaller than those required for the adiabatic 
approximation [see Eq. (17)], novel features appear in 



the exact diagonalization approach such as the possibil- 
ity to simulate a spin-orbit coupling between two Bosc- 
Einstein condensates in different dressed atomic states, 



experimentally observed by Lin et al. (2011). 



The discussion of Sec. II. El about the observable vortex 
pattern applies directly to the present scheme. The num- 
ber of vortex rows that can fit along the y axis is again 
given by ~ 4\/Pv ~ V^A- The nim t of large vortex 
arrays can thus in principle be reached with this config- 
uration, provided one uses a relatively small gradient of 



detuning. Lin et al. (2009a I generated about 10 vortices 
in a rubidium Bose-Einstein condensate. These vortices 
did not order into a regular lattice, presumably because 
of the heating of the cloud and atom losses (1/e lifetime 
of 1.4 s) caused by the residual photon scattering. 

One could think that increasing the detuning of the 
laser coupling would solve the photon scattering prob- 
lem. Unfortunately the situation is not so favorable, at 
least for alkali-metal species. Indeed, by contrast to the 
scalar light shift that scales as U sca \, oc I/A where / 
is the light intensity, the two-photon Raman coupling 
k is proportional to the vector part of the light shift 



U vec . oc /Aps/A 2 . Here Aps denotes the fine structure 
splitting, which is only a few percents of wa for alkalis 
(2% for rubidium atoms, 5% for cesium atoms) and the 
above scaling for k is valid for Aps ^ A. Therefore in the 
limit of large detunings the Raman coupling k decreases 
with A as fast as the photon scattering rate cx IT/ A 2 . 
In other words, the usual trick that consists in increasing 
both A and I to enhanced the role of the scalar light 
shift with respect to spontaneous emission processes is 
not applicable in this context. 



III. NON-ABELIAN GAUGE POTENTIALS 

The idea of non-Abeli an geometric gauge pote ntials 
goes back to the work by Wilczek and Zee (1984), who 



considered the generalization of the adiabatic theorem to 
the case where the Hamiltonian of the system of interest 
possesses a group of eigenstates that remain degenerate 
(or quasi-degenerate) and well-isolated from other levels 
in the course of the time evolution. This analysis was fol- 
lowed by applications to many areas including inter alia 



molecular and condensed matter physics (Bohm et al. 



2003 Xiao et al. 20101. In particular, the possibility to 
generate non-Abelian magnetic monopoles was demon- 
strated in the rotational dynamics of diatomic molecules 



(Bohm et al. 


1992 


Moody et al. 


1986 


Zygelman 


L990) 


and in the nuclear quadrupole resonance ( 


Zee 


198E 


>). 



In this section we are interested in the non-Abelian 
dynamics of cold atoms in light fields, which can emerge 
when a (TV + l)-state atomic system with N > 3 (see 
Fig. [6]) is excited by a suitable configuration of laser 
beams. The first study in this direction was performed 
in the context of Stimulated Raman Adiabatic Passage 
(STIRAP) by |Unanyan et al.\ ( [19981 11999[ ), who consid- 
ered the behaviour of a four-state atomic system (tripod 
configuration) when it is driven by three successive laser 
pulses . Subs equently Osterloh et al. (2005) and Ruseckas 
et al.\ ( |2005 ) transposed the concept of gauge potentials 
to the case of a continuous atom-laser excitation relevant 
for the present review and identified situations where the 
non-Abelian gauge potentials emerge. 

We start this section by providing a general formula- 
tion of the adiabatic motion of atoms when some internal 
atomic states remain degenerate in the presence of the 
atom-light coupling. We show how non-Abelian gauge 
potentials appear and discuss the structure of these po- 
tentials for some typical laser configurations. We also 
present some physical phenomena that are associated 
with non-Abelian gauge potentials, such as the genera- 
tion of magnetic monopoles, Rashba spin-orbit coupling, 
and non-Abelian Aharonov-Bohm effect. 



A. Emergence of non-Abelian gauge potentials 

In this section we consider atoms with N + 1 inter- 
nal levels, which are coupled to a light field. We assume 



11 



A Emergence of non-Abelian gauge potentials 



III NON-ABELIAN GAUGE POTENTIALS 



that we can cast the atom-laser interaction at any given 
point r as a time-independent (N + 1) x (N + 1) ma- 
trix U (r) using the rotating wave approximation. Our 
starting point is then still the full atomic Hamiltonian of 
Eq. ([!]). For a fixed position r, U(r) can be diagonal- 
ized to give a set of N + 1 dressed states \Xn(r)), with 
eigenvalues e n (r) (n — 1, . . . ,N + 1). The key feature 
that leads to non-Abelian gauge potentials is that some 
dressed states form a degenerate (or nearly degenerate) 
manifold, at any point r in space. More specifically we 
shall assume that the first q atomic dressed states form 
a degenerate subspace £ q , and these levels are well sepa- 
rated from the remaining ones. 

The full quantum state of the atom describing both 
internal and motional degrees of freedom can be written 
I*) = Eri^W \Xn(r)) where each tp n is the wave- 
function for the centre of mass motion of the atom in 
the internal state \xn)- We are interested here in the 
dynamics of the atom when it is initially prepared in the 
subspace £ q . Neglecting transitions to states outside £ q , 
we can project the full Schrodinger equation onto £ q and 
we arrive at a closed Schrodinger equation for the reduced 
column vector ^ = . . . , ipq) 1 



lh -dt = 



2M 



+ (V + e)l q + W 



(38) 



where l q is the identity matrix in £ q and e(r) is a diag- 
onal matrix of eigenenergies £ n {r) (n — 1, .. .,<?). This 
equation is very reminiscent of Eq. ([6]); However A and 
W are now q x q matrices with the elements 



^ N+l 

Wn.m = — V! A n , ■ At. 



2M 



(39) 
(40) 



2=9+1 



with n, m € (1, . . . , q). The effective vector potential A 
is called the Mead-Berry connection. 

The effective magnetic field (or curvature) B associ- 
ated to A is: 



1 



Uki Fui, Fh, = dhA, — d,Ah 



■[A k ,A t ]. (41) 



Note that the term ^Eiki[Ak, A{\ = (A X A)i does not 
vanish in general, because the vector components of A 
do not necessarily commute. Therefore the magnetic field 
B can be non-zero even if the vector potential A is uni- 
form in space. This property is specific of non-Abelian 
dynamics and occurs only if q > 2, whereas for q = 1 
we recover simply the Abelian dynamics discussed in the 
two preceding sections. 



B. The multipod scheme 

A generic way to obtain a degenerate subspace in atom- 
laser interaction is to realize the situation depicted in 




FIG. 6 Multipod configuration. An atomic state |e) is cou- 
pled to N different atomic states \gj) (j = 1,...,N) by N 
resonant laser fields. 



Fig. [6j where one single level labelled |e) is coupled to 
N levels labelled |g 7 -) (j = 1, . . . , N), with complex Rabi 



frequencies Kj (Juzeliunas et al. 
at the end of this subsection he 



implemented in practice for N 
(quadrupod). The atom-light coupling operator is 



2010J). We will discuss 
n this scheme can be 
3 (tripod) and N = 4 



y, ft/Cj(r) 

^ 2 

.7=1 



eXffjI+H.c. 



which can be conveniently rewritten as 



U = 



hn{r) 



(\e)(B(r)\ + \B(r))(e\) 



(42) 



(43) 



where \B) = YljLi K j\9j)/ K i s the bright (coupled) state 
generalizing ( 26 ) and k is the total Rabi frequency, k 2 — 

The diagonalization of U is straightforward. First 
the coupling between \B) and |e) gives rise to the two 
eigenstates |±) = (|e) ± \B)) /V2 with energies ±Hk/2. 
Then the remaining orthogonal subspace of dimension 
N — 1 corresponds to degenerate dark states that are 
all eigenstates of U with energy e = 0. This provides 
the degenerate subspace £ q with q = N — 1 introduced 



in Sec. III.A[ which is required for the emergence of 
non-Abelian dynamics. In the following we denote \D n ) 
(n = 1, . . . , N—l) a normalized orthogonal basis of £n-i- 
Let us briefly discuss how such a multipod scheme 
can be implemented in quantum optics for the two cases 
N = 3 and N = 4. The scheme of Fig. [6] for N = 3 is re- 
alized in a straightforward way by considering an atomic 
transition between a ground electronic state with angu- 
lar momentum J g = 1, and an excited electronic state 
with angular momentum J e = 0. Such a transition oc- 
curs for alkali-metal species such as 23 Na or 87 Rb, as well 
as for helium atoms prepared in metastable electronic 
spin triplet states ( 3 <Si — 3 Po transition). The scheme 
N = 4 is a bit more subtle to achieve and we briefly 
outline the proposal detailed by Juzeliunas et al. (2010). 



The idea is to use the two hyperfine ground levels of 
an alkali-metal atom like 87 Rb, with angular momenta 
equal to F = 1 and F = 2, respectively (see Fig. 
For the state |e) we choose one particular ground state 
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(b) 



F = 2 



F = 1 



FIG. 7 Implementation of the multipod scheme for TV = 4 
for alkali-metal atoms with two hyperfine levels of angular 
momentum F = 1 and F = 2. The couplings involved 
in this scheme correspond to stimulated Raman transitions 
between hyperfine states of the ground levels. We choose 
|e) = \F = = 0). A small magnetic field lifts the 

degeneracy within the levels F = 1 and F — 2 and al- 
lows one to address selectively the various Raman transi- 
tions, using pairs of laser beams with properly chosen po- 
larization and frequencies, (a) Two pairs of laser beams in- 
duce the transitions |e) — > \gi) = \F = 2, mf = 1} and 
e) — > 1^3} = \F — 2,niF = — 1). (b) A pair of laser beams 
induce the transitions |e) — > \g2) = \F = l,m,F = 1} and 
\e)->\g4) = \F = l,m F = -l). 



\e) = \F = l,m = 0). The four other states play- 
ing the role of the states \gj) are the Zeeman sublevels 
\F = 1, mi? = ±1) and \F = 2,mp = ±1). The couplings 
between |e) and the levels \gj) are induced by resonant 
two-photon Raman transitions like in Sec. |II.D| The de- 
coherence and heating due to photon scattering can be 
minimized in this scheme by chosing a large single pho- 
ton detuning (that is, of the order of the fine structure 
splitting Aps), as explained in Sec. II. D 



in which case we can choose the following expression for 
the dark states in the basis {|<?i), I32), |<?3)} 



\Di) = 



(cos a cos j3 e^ 31 ' 
cos a sin/3e 4 * 32 
— sin a 



The vector potential then reads 



(45) 



sin 2 (3 V<; 



'13 J 



with iftij - 

An = h (cos 2 (3 V 

A 22 = h cos z a (cos^ P V0i3 + sin^ /3 V 



1 



/icosa ( -sin(2/3) - ) (46) 



'23 



To generate a magnetic monopole, we consider the 
laser setup formed with two of the laser beams propa- 
gating along the z axis, and the third one along the x 
axis. More precisely the first two beams are prepared in 
the Laguerre-Gauss modes with orbital angular momen- 
tum I = ±1, and the third beam is prepared in the first 
order Hermite-Gauss mode so that 



M,2 = Ko|e i ^), 



K3 



n -e 



(47) 



Here p is the distance from the z axis and (p the azimuthal 
angle around this axis. We suppose that the beam waists 
are large compared to the other physical scales of the 
problem and we therefore omit the Gaussian profiles in 
Eq. (47). The calculation of the vector potential using 



the formalism of Sec. |III.B| gives 



A = 



tantf^ 



hk 



(e z - e x ) [(1 + cos 2 i?)l + (1 - cos 2 d)a z ] , 



C. Generating a magnetic monopole 



where r, z9 and ip are the spherical coordinates, &j are the 
Pauli matrices and 1 is the 2x2 unity matrix. The first 
term in the vector potential proportional to a x represents 
the field of a magnetic monopole of unit strength at the 
origin (r = 0): 



One of the interesting properties of non-Abelian 
gauge potentials is the possibility to generate magnetic 



monopoles, as first pointed out by Moody et al. (1986). 



cment 



In this section we present a possible way to imp 
such a monopole using atom-light interaction. We con- 
centrate on the physical aspects of the problem and we 
refere the reader to the original publication by |Ruseckas 
et al.\ pOOBl for technical details. 



We will consider here the tripod configuration obtained 



using N = 3 ground levels in the formalism of Sec. III.B 



It is convenient to parametrize the Rabi frequencies Kj 
with angle and phase variables according to 



Ki = k sin a cos j3 e 1 
K2 = n sin a sin /3 e 1 ^ 2 



(44) 



«3 



k cos a e 



*f>3 



B = 



6 r (J . r 



(48) 



The dots indicate non-monopole field contributions 1 pro- 
portional to the Pauli matrices and to the unity matrix. 



Note also that the total intensity of the laser fields ( 47 ) 



vanishes at the origin representing a singular point. Thus 
one should apply an additional potential that will ex- 
pel the atoms from this region, in order to avoid non- 
adiabatic transitions in the vicinity of r = 0. 



A pure monopole would emerge if the Rabi frequencies Kj could 
be taken proportional to the corresponding three Cartesian co- 
ordinates ki = Ax, K2 = Ay and K3 = Az. This is however not 
possible in practice, as this spatial dependence is not consistent 
with the Maxwell equations for the driving laser field. 
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This configuration was further analyzed by |Pietila and 
Mottonen (2009), who studied the behavior of an inter- 



acting Bose Einstein condensate in this monopole con- 
figuration. They showed numerically that the existence 
of the monopole gives rise to a pseudospin texture with 
a topological charge that cancels the monopole charge. 



D. Generating a spin-orbit coupling 

Non-Abelian light-induced gauge potentials can be 
used for generating a spin-orbit coupling for cold atoms, 
simulating the one appearing for electrons in condensed 
matter. The electron's spin degree of freedom plays a 
key role in the emerging area of semiconductor spintron- 

A 



Fert 


2008 


Griinberg 2008 Zutic et al. 


2004) 



first scheme for a semiconductor device is the spin field- 
effect Datta- Das transistor (DDT). It was proposed 20 
years ago ( Datta and Das| |1990[ ) and implemented re- 
cently ( Koo et al. 2009 ) . Atomic and polaritonic analogs 



of the electron spin transistor have also been suggested 
(Johne et al. 2010 Vaishnav et al. 20081. An impor- 



tant ingredient of the DDT is the spin-orbit coupling of 
the Rashba (]Rashba[|1960||Wmkler||2003|) or Dresselhaus 
(Dresselhaus 1955 Schliemann et al. |2003 ) types. This 
Rashba-Dresselhaus (RD) coupling scheme is described 
by a vector potential which can be made proportional to 
the spin- 1/ 2 operator of a particle within a plane ( Cserti 
and David 



2006) 



We explain in this section how an effective spin-orbit 
coupling can be generated for atoms for the cases of an 
effective spin s = 1/2 and s = 1. We start with the gen- 
eral formalism of a iV-pod level scheme in the case where 
the coupling lasers are plane waves of equal amplitude 
propagating in the xy plane. More precisely we assume 
that the wave-vectors kj form a regular polygon 



k(- 



■e x cos CYj 



e v sm a 



), aj = 27rj/N. (49) 



We set Kj = k^ ) e xki ' r /\N and use the set of orthogonal 



normalized dark states 

N 



,N-1. 



\D n ) = -^J2\ 9j )e i ^ n - ik r-, n = l, 

(50) 

Substituting \D n ) into Eqs. (39)-(40), one arrives at con- 



stant (yet non-Abelian) scalar and vector potentials 
H 2 k 2 



AM 
hk 

~V2 



(5 m ,l8n,l + dm,N-l8n,N-l) 7 (51) 



n,m±l 



(52) 



with e± = (e x ± ie y )/\/2. For a constant exter- 
nal potential V, a basis of stationary solutions of 
the Schrodinger equation (38) are the plane waves 
&K(r,t) = $ K - e i ( JC ' r-nK " t ) with the amplitude <S> K 



obeying the eigenvalue equation Hk&k = ^k^k 
where Hk is the i^-dependent (N — 1) x (N — 1) matrix: 



H K = 



(HK1 



2M 



W + Vl 



(53) 



In the tripod setup (TV = 3) the wave vectors kj form 
an equilateral triangle and W, A and Hk are 2x2 ma- 
trices. The scalar potential is proportional to the unit 
matrix, W = h 2 k 2 /(AM) 1, whereas the vector poten- 
tial A = —kS± is proportional to the spin 1/2 oper- 
ator S± = S x e x + S y e y in the xy plane (S = her/2, 
where the <x are the Pauli matrices). This provides 
a spin-orbit coupling of the RD type, characterized by 
two dispersion branches ft K = h(K ± k/2) 2 /2M for 
V = —h 2 k 2 /(AM). A number of other arrangements of 
laser beams have been considered to produce the same 



RD spin-orbit coupling (Jacob et al. 


2007 


Juzeliunas 


et al.\ 2008b Larson and Sjoqvist, 2009 Stanescu et al.\ 


2007 


Vaishnav and Clark 


2008| Zhang et al. 


2010 


)• 



In the tetrapod setup (N = 4) the choice of Eq. (|49j) 
corresponds to two orthogonal pairs of counterpropagat- 
ing laser fields. The vector potential A = —kJ±/\/2 is 
proportional to the projection of a spin 1 operator in the 



xy plane J± = J x e x 



Jy e v> 



whereas the scalar poten- 



tial is proportional to the squared z component of the 
spin operator, W = J 2 k 2 /(AM). The eigenfrequencies 
are now 

2M { 



nn p K = 



K z 



2Kk/3+k 2 )+V , /3 = 0,±1, (54) 



with K = \K\. For j3 = ±1 the dispersion curves are 
analogous to those of the spin- 1/2 RD model. The ad- 
ditional dispersion curve with /3 = corresponds to a 
parabola centered at K = 0. 

A spectacular consequence of spin-orbit RD coupling is 
the Zitterbewegung, a phenonenon which was analyzed 
for cold atoms (|Larson et al.\ |2010| |Merkl et al.\ |2008 



Song and Foreman| |2009 Vaishnav and Clark 



electrons i n solids^Cserti and David 2006 Rusin and 



Zawadzki 2009 



2008) 



Schliemann et al.\ 2006 ) and trap ped 
ions ( |Lamata et cd. 2007 Rusin and Zawadzki 20101. It 
was recently observed for the latter systems ( |Gerritsma 



et al. 2010). Another manifestation of the RD cou- 
pling is the negative refraction and reflection that occurs 
when a matter wave is incident on a potential step. This 
problem was investigated for spin- 1/2 atoms (Juzeliunas 



and Winkler 



et al. 2008a) and electrons (Dargys 2010 Teodorescu 



2009) 



For small incident wavenumbers, 
K <C k, the transmission probability is close to unity 
for normal incidence. This nearly complete transmission 
is a manifestation of the Klein paradox appearing also 
for electron tunneling in graphene (|Katsnelson et ~aL 



2006 Neto et al. 2009[ ). The transmission probability 



decreases when the angle of incidence increases, and the 
transmitted matter wave experiences negative refraction, 



again similar to electrons in graphene (Cheianov et al. 
2007). Particles with a spin larger than 1/2 have ad- 



ditional internal degrees of freedom, which modifies the 
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continuity conditions at the potential step. For the neg- 
ative refraction phenomenon at non-normal incidence, 



Juzelhmas et al. (2010) have shown that the amplitude 



of the refracted beam is significantly increased in com- 
parison with the spin 1/2 case. 

Spin-orbit coupling with its multicomponent disper- 
sion curves also offers interesting scenarios in the pres- 



ence ofcollisixjnal interactions ( Stanescu et al 



2007 Wang et al. 2010 Yip 2010). In two dimen 



2008 



sions with a symmetric RD gauge potential of the form 
hk(a x e x +a y e y ) , the dispersion curve has a minimum at a 
constant non-zero radius in momentum space. This cor- 
responds to a massive degeneracy of the single-particle 
ground level, since all plane waves with momentum K 
such that \K\ — k are possible ground states. For a Bose 
gas, this precludes in principle the formation of a conden- 
sate, and atomic interactions will lead to the formation of 
a strongly correlated ground state. However in the pres- 
ence of an asymmetry in the RD coupling, the massive 
degeneracy is lifted and the minimum of the dispersion 
occurs for two opposite values of the atomic momenta 



(Stanescu et al. 2008 20071. Weakly-interacting bosons 



with such RD coupling will then behave essentially as a 
two-component system at low temperature. 

We conclude by signalling a possible drawback of the 
tripod scheme: The two dark state states forming the ef- 
fective spin-1/2 system are not the lowest single-particle 
energy states, hence collisional deexcitation can trans- 
fer the atoms out of the dark state manifold down to 
the ground dressed state. To overcome this difficulty, a 
scheme involving N ground or metastable internal states 
cyclically coupled by laser fields was recently proposed 
(Campbell et al. 2011b. By properly setting the direc- 



tion and phases of the laser fields, a pair of degenerate 
pscudospin states with lowest energy emerge. The states 
are subjected to the RD coupling and are immune to the 
collisional decay. 



E. Non-Abelian Aharonov-Bohm effect 

The spin-orbit coupling that we presented in the 
preceding subsection can also be the source of quasi- 
relativistic dynamics, where the atomic centre of mass 
motion is governed by an effective Dirac equation in the 
limit of small momenta and strong gauge potentials. The 
quasi-relativistic dynamics of the atomic centre of mass 
is a remarkable effect which connects phenomena from 
high energy physics, graphene and photonic crystals with 
atomic physics and quantum optics. It relies on the two- 
or three-component nature of the effective atomic inter- 
nal structure, but it is not explicitly linked to the non- 
Abelian nature of the coupling. 

To observe a non-Abelian effect one should study spin 
dynamics. In our case it will be a pseudo spin, i.e. an 
effective multi-component system created by the interac- 
tion between incident laser beams and the atoms. Pre- 
viously we have shown that the gauge potential affecting 



the pseudo spin can be for instance proportional to the 
Pauli spin matrices in the x and j/-direction such that 
the matrices corresponding to the two directions do not 
commute. This is indeed also the requirement for having 
a non-Abelian system. A simple illustration of a non- 
Abelian situation is the so called non-Abelian Aharonov- 



Bohm experiment (Jacob et al. 2007). Here we envisage 



a particle with an internal structure, the pseudo spin, 
which is allowed to move along two different paths from 
A to B along straight lines in x and y-direction. Let 
us suppose the path consists of two equidistant trips of 
length L along the x and y directions. The question is, 
if we start at A, with a given orientation of the pseudo 
spin, what is the spin at the final point B? If the spin 
is subject to a non-Abelian field, the result will depend 
on which path the particle took. This is easily seen by 
neglecting any external dynamics and only considering 
the effects from the gauge potential, which is of the form 
A = A, x e, x + A y e y . The final state based on the path 
going first along x-direction and then y-direction will be 
given by exp[zA y L] exp^A^Lj^in where 'J'in is the initial 
pseudo spin. If we take the path along the y-direction 
first and then in the cc-direction, we will have the final 
state exp[iA x L] ex-p^AyL^^. These two states are not 
necessarily the same, because A x and A y do not com- 
mute if the system is non-Abelian. For instance, when 
A x oc a y and A y oc a x , these operations are equivalent 
to spin rotations around x and y, which do not yield the 
same final state in general. Interestingly this situation is 
similar to the scattering of protons onto a non-Abelian 
flux line where the protons are anticipated to be con- 



verted into neutrons (Horvathy 19861 



IV. GAUGE POTENTIALS IN OPTICAL LATTICES 

Optical lattices have recently emerged as a major tool 
for th e fie ld of quantum gas es [see Lewenstein et al. 
d2007[) and|Bloch et aZ.|(|2008[) for recent reviews]. They 



correspond to a periodic array of trapping sites connected 
to each other by quantum tunneling. Such a potential 
landscape is created using the interference pattern of sev- 
eral off-resonant lasers. In addition to the obvious anal- 
ogy with the effective periodic potential exerted by the 
ionic matrix on electrons in a real solid, optical lattices 
allow one to bring quantum gases into the strongly corre- 
lated regime, where interactions dominate the behavior of 
the system. The entrance in the correlated regime is ob- 
tained by increasing the lattice depth and it follows from 
two effects, the increase of the on-site interaction energy 
due to the stronger confinement of the atoms near the lat- 
tice sites, and the reduction of the tunneling probability 
from site to site which decreases the kinetic energy. 

In this section we first review some basic elements on 
the band structure in a periodic potential (§ |IV.A I and 
the effec ts of a magnetic field on the single particle spec- 
trum (§ IV. B). We then discuss the concept of laser as- 



sisted tunnelling, which allows one to control the phase of 
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the tunneling matrix elements and realize artificial gauge 
potentials (§ 



position r„. m and the initial Hamiltonian P 2 /(2M) + V\ a 



IV. C and IV.D ). Very recently another class can be replaced by 



of lattices with orbital magnetism, the so-called flux lat- 
tices, has been proposed by Cooper ( |2011 1. We briefly 
relate this proposal as well as those based on laser as- 
sisted tunnelling to the dressed state approach discussed 
in the previous section in § |IV.E| Finally we give a short 
account on the possibility to generate non-Abelian fields 
in lattices in § |IV.F| 

Let us mention briefly some alternative methods that 
do not use laser coupling to realize effective magnetic 
fields for cold atoms in optical lattices. One is simply 
to rotate the lattice, which was realized experimentally 

( |2010} for lat- 



H 



TB — — J ®n,m®n',m' + h.C, 



(57) 



where a nm is an annihilation operator for a particle in 



the state 



The sum runs over nearest neighbors 



Tung et al. ( 2006 ) and Williams et al 



tices with large lattice spacing, i.e. outside the Hubbard 
regime. Another method relies on a temporal modulation 
of the la ttice potential wit h x—y ( S0rensen et al 2005|) or 



x" — y ( Lim et al. |2008 1 symmetries (see also |Kolovsky| 
(2011)). One can also control the sign of the tunnelling 



matrix element by modulating of the lattice potential 
Eckardt et al.\ (|2005|, |Lignier et~aT\ (|2007[)], which can 



provide artificia 



magnetism for some specific lattice ge- 



ometries. Finally on can obtain the desired gauge poten- 



(n.n.) only, and J is the tunelling energy character- 
izing the hopping between neighboring sites. Hopping 
to more distant sites is neglected in this approximation. 
The single-particle dispersion is in the TB approxima- 
tion e(q) — —2J{cos(q x d) + cos(q y d)], corresponding to 
an energy width of 8 J for the band. 



B. Harper equation and Hofstadter butterfly 

Suppose now that the square lattice considered above 
is placed in a uniform magnetic field B = B e z , asso- 
ciated to the vector potential A — (—By, 0,0) (Landau 
gauge). For a particle carrying a charge e, the Hamilto- 
nian (57) can be written 



tial through interactions with a bath of atoms (Klein and 



Jaksch 2009) 



u 'n±l.m u "n,,m 



A. Reminder on band structure 

Consider a particle moving in the two-dimensional 
square lattice potential with period d and depth Vo 

V iat (x, y) = Vo [sin 2 (irx/d) + sin 2 (ny/d)] . (55) 

The energy eigenstates are the Bloch waves Vyq( r ) 



n,m,± 

(58) 



n,m,± 



with 



A ■ ds, 



(59) 



where the integral is taken along a straight line joining 
neighboring sites. The appearance of the phase factors 



^(r), where n = 0, 1, . . . denotes the band index e±i *„, m in Eq pi C£m be understood in terms of the 



and the quasi momentum q takes values in the first Bril 
louin Zone (BZ) | — n/d, ir/d\x \— ir/d, ir/d\ [see for exam- 



ple Ashcroft and Mermin (1976)]. The function u Vtq is 



a periodic function of r, with period d along both direc- 
tions of space. The energies e n (q) associated with the 
Bloch waves form bands when q is varied inside the first 
Brillouin zone. 

We assume in the following that the lattice is in the 
tight-binding (TB) regime, where the energy width of 
each band is much smaller than the gaps between two 
consecutive bands. More specifically, we consider par- 
ticles confined to the lowest Bloch band 77 = as it is 
usually the target regime for experiments on correlated 
quantum gases, and we subsequently drop the band in- 
dex. In this limit, it is useful to transform to the orthog- 
onal basis of Wannier functions 



Aharonov-Bohm phase accumulated along a straight line 
joining two adjacent sites of the lattice. With the choice 
of the Landau gauge, these phase factors are 



n+l,m 



2nam, 



(60) 



where a = $ = Bd 2 is the magnetic flux through 

a unit cell and $0 = h/e is the flux quantum. A different 
gauge choice would lead to different phase factors 4> n ,rm 
but the phase accumulated on a closed circuit around an 
elementary cell 



7 



eBd 2 /h=2na 



(61) 



t {r)=Af 



BZ 



%{r) d 2 q (56) 



where Af is a normalization factor and r n , m = d{ne x + 
me y ), with n,m being integers. In the TB regime, the 
Wannier function w n m is localized near the lattice site at 



is gauge invariant (the symbol □ indicates that the sum 
runs over an elementary cell of the lattice). 

The reduction of the full Schrodingcr equation to 
Eq. ( 58 ) involves another approximation, the so-called 
Peierls's substitution [see Luttinger ( |1951[ ) for the dis- 
cussion of this approximation in the TB limit] . It is valid 
if the Landau energy huj c , where lo c = B/M is the cy- 
clotron frequency, remains much smaller than the energy 
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1/2- 



-L 



FIG. 8 Hofstadter butterfly: single particle spectrum of the 
Hamiltonian of Eq. ( 58 1 in units of the tunneling amplitude 
J, for a ■ 




varying between and 1. 



gap between the bands n = and n = 1. Outside of 
the TB regime, its validity has to be carefully examined 
(Nenciu 1991). In the context of our Colloquium, these 



considerations are not relevant as we will discuss methods 



to simulate the Hamiltonian ( 58 1 directly. 



The single particle eigenstates of the Hamiltonian ( 58 ) 
can be searched in the form ^ ran C ra e I,I "' i oj lim |vac). 
The corresponding eigenvalue equation for the coeffi- 
cients C rn is known as the Harper equation ( |Harper[ 
1955 ) . Its solutions have been extensively studied in con- 



nection with the quantum Hall effect (Kohmoto 1989 



in Fig. 8 



Thouless et al. 1982). The energy spectrum, represented 



is known as the Hofstadter butterfly (Hofs- 
tadter 1976), and it exhibits a remarkable self-similarity. 



This spectrum has a much more complex structure than 
its counterpart for a continuous system, the latter be- 
ing composed of uniformly spaced, infinitely degenerate 
Landau levels. This structure can be qualitatively un- 
derstood on a relatively simple ground. For a ratio- 
nal value of a, that is a = p/q with p, q integers, the 
Hamiltonian is still periodic but with a larger period qd 
along the y axis. One can then divide the system into 
"macro-cells" of size qd and look for new eigenstates [the 
so-called magnetic Block functions (Blount 1962J)] ap- 
propriate for the "macro-lattice". Because the macro- 
cell contains q original lattice sites, each energy level is 
q— fold degenerate and the fundamental band splits into 
q subbands separated by small energy gaps (with the ex- 
ception of a = 1/2, where the two bands touch at the 
edges of the reduced Brillouin zone). The Landau lev- 
els structure is recovered in the weak flux limit where 
a = p/<7 <C 1. We note finally the special case of a = 1/2, 
corresponding to a real tunneling matrix element alter- 
nating in sign from one column to the next. In this 
case, the excitation spectrum features two 'Dirac points' 
for (k x ,k y ) = (±ir /d x ,ir /2d v ) around which the disper- 
sion relation is linear (Hatsugai and Kohmoto 1990 Hou 
et al\ |2009| |Lim et al.[ |2008[ ). The behavior of an ul 
tracold Fermi gas in such a situatio n is expected to b e 
similar to that observed in graphene ( Neto et al. 2009 1 . 



y/d y ^ 
m + 1 



in 



m — 1 



-Q- 



i 



2ira 



-O 



-Q- 



t 



p— - x/d x 



n-1/2 n n + 1/2 n + 1 

FIG. 9 (Color online) Sketch of the state-dependent 2D lat- 
tice potential. Grey (black) dots mark lattice sites for atoms 
in the ground g (excited e) state, respectively. We highlight a 
trajectory around an elementary cell of the coupled lattice in 
which the particle acquires the Aharonov-Bohm phase 2na. 
The g — > g and e — > e transitions occur thanks to standard 
tunneling. The g — > e and e — > g transitions correspond 
to resonant laser assisted tunneling. In a g — > e transition 
the atom absorbs a photon and jumps from a state localized 
around x — nd x to a state localized around x = (n ± l/2)d x 
(and vice versa for a e — > g transition). 



C. Simulating a magnetic flux through each lattice cell 



We now discuss how the Hamiltonian (58) can be re- 
alized for cold atoms in an optical lattice. We proceed in 
two steps. First we introduce the notion of laser assisted 
tunneling and show how it can be used to obtain a non- 
zero flux through each lattice cell. However laser assisted 
tunneling in its simplest version provides a flux that al- 
ternates in sign between neighboring cells, and thus does 
not simulate a uniform magnetic field. In a second step 
(Sec. IV.D), we will show how one can rectify this mag- 



netic flux and reach the Hamiltonian ( 58 1 



The notion of laser assisted tunneling has been intro- 
duced in this context by Ruostekoski et al. (20021 and 
Jaksch and Zoller (2003). The first ingredient is to de- 



sign a state- dependent lattice. Consider atoms with two 
internal states g and e, trapped in spatially separated 
sublattices (see Fig. [9| . We focus on the situation where 
each sublattice has rectangular symmetry with lattice 
spacings d x and d y . The e sublattice is deduced from 
the g sublattice by a translation d x /2 along the x axis. 
We label 



2n,m 
2n+l,m 



(n+ 1/2) + mdyCy . 



(62) 
(63) 



the positions of the trapping sites on each sublattice, with 



the corresponding Wannier functions w^nm an< ^ 



(e) 
' J 2n+l,rn- 



Tunneling energies along x and y within a given sublat- 
tice (g or e) are denoted by J x and J y , respectively. 

In practice, there are two options to generate such a 
state-dependent lattice. The first one, suitable for alkali 
atoms, selects g and e as two Zeeman or hyperfine states 
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in the electronic ground state manifold (electronically ex- 
cited states are too short lived for this application). In 
this case, the leading terms in the optical lattice potential 
is the sum of scalar and vector terms, the latter behav- 



ing as an effective magnetic field ( Cohen- Tannoudji and 



Dupont-Roc 1972). The states g and e are chosen such 
that they have opposite magnetic moments, and the de- 
tuning from resonance of the light beams forming the x 
lattice is adjusted to that the polarizability of g and e are 



opposite (Mandel et al. 2003). This provides the desired 
state dependent lattice potential along the x axis. Note 
however that the required detuning lies relatively close 
to the resonance lines, so that heating due to photon 
scattering can cause serious practical problems. 

The second option is to use atomic species already con- 
sidered in Sec. |LB with a long-lived excited state ( Gerbier 



and Dalibard 2010 Yi et al. 20081. Practical examples 
are alkaline-earth atoms or Ytterbium atoms, where the 
3 Po internal state has a typical lifetime over ten seconds. 
One then chooses the electronic (spin singlet) ground 
state for g and the 3 Jo excited state for e. For such 
atoms, both states couple to different electronic states, 
and it is possible to find a magic wavelength (used for 
the y lattice) such that the polarizabilites are identical 
for both states, and an anti-magic wavelength (used for 
the x lattice) such that they are opposite. Both are far 
detuned from any resonance, so that heating by sponta- 
neous emission is not an issue in this case. 

In a state-dependent lattice, tunneling from a given site 
on the g sublattice to a neighboring one on the e sublat- 
tice can be driven by a resonant light field that couples 
g and e. We assume this coupling laser is a plane run- 
ning wave with wavevector k. The corresponding com- 
plex tunneling matrix element is 



a r (s) 

Hi ' 2n,m 



-> e, *4 e > : J cS = Sk<°> O e ik -<™, (64) 



where characterizes the strength of the atom-laser 
coupling and the dimensionless number 



O 



(«41i, m (r)) «fijr)e*'A (65) 



is the overlap integral between ne ighboring Wannier func- 



tions (Jaksch and Zollcr 

(e) 

2ra+l,m 



20031) (note that w^ m and 



w 



are not orthogonal since they belong to different 



sublattices). Eq. (64) displays the essential phase factor 



that allows to reproduce the Aharonov-Bohm phase (59 ). 
By choosing a laser propagating in the y — z plane, and 
under the assumption that state independent tunneling 
along x is negligible (J x <C \ J c g\), one obtains 



H = 



\Je S \ £ 



i<t>, 



n,m,± 



^L+l,m±1^2n+l,mt 



(66) 



n,m,± 



where a\ (resp. b 



t 

2n+l,m 



) creates an atom in the inter- 



nal state g (resp. e) in the state w-fnm ( res P- w 2n+i 



„( e ) 



Here e l ^ m is the complex phase of J c s, characterizing the 
phase acquired by an atom when it tunnels from a site 
of the sublattice g to a site of the sublattice e: 



a r (s) 

ill ' 2n.r 



— > e, r 



(e) 

2n± X,m 



27ram 



(67) 



where a = kydy/(2Tr) can be adjusted between and 1 by 
changing the angle of the wave vector k of the coupling 
laser with the y axis. 



Although the Hamiltonian ( 66 ) contains complex hop- 
ping amplitudes, it does not yet coincide with the Hamil- 
tonian ( 58 1 that we wish to simulate. In the initial Hamil- 



tonian, the phase picked up across each ic-link has the 
same sign for a given link direction in real space: 2itam 
in the +x direction (and thus — 2i:am in the —x direction 
since the Hamiltonian is Hermitian). Here, the phase is 
tied to the direction in internal space: 2itotm for transi- 
tions to g to e, irrespective of the fact that they occur in 
the +x or —x directions. 



In total, we achieve with ( 66 ) a situation where the flux 



through a cell changes sign from one cell to the next one 
along the x axis. We anticipated this difficulty in Sec |I.E| 
where we noted that realizing artificial gauge potentials 
with variations on the scale of the optical wavelength 
Xl would naturally result in an oscillating (or staggered) 
effective magnetic field. Although this configuration can 



lead to interesting situations 1 


Hou et al. 


2009 


Moller 


and Cooper 


2010 


Wang and Gong 


2006), it does not 



realize our goal of producing a uniform magnetic field, 
except for the specific case a — 1/2, since phase changes 
of ±7r are equivalent. 



D. Rectification of the magnetic field in the lattice 

The need for rectification of the magnetic field that 
we found in the preceding subsection can be formulated 
on more general grounds. We are looking for an effec- 
tive Hamiltonian that must break time-reversal symme- 
try. For our t wo-level system, the time reversal operator 



is T = K. S K. C (Messiah 
is such that lC s \g) = 



, where IC S = eJ 



,/2 = 



1961a) 

|e), JC s \e) = —\g), and JC C is the 
conjugation operator. Now, the action of T on the con- 
figuration studied in the previous section results in the 
same configuration, except for a mere translation by an 
amount d x /2 along the x axis. Hence, in order to pro- 
duce a uniform magnetic field, one must add another 
ingredient that is not symmetric under the "sublattice 
translation" operation. 



|Jaksch and Zoller (2003) proposed to realize this by 
adding a linear potential gradient along x, V^ ra d(n) = (n, 
so that the laser frequency for the (<?; 2n, m) — > (e; 2n + 
1, m) transition becomes shifted by +£ while that for the 
(<?; 2n + 2,m) — > (e; 2n + l,m) transition is shifted by 
— £• By choosing ( ^> \J e e\, the two transition frequen- 
cies become non-degenerate and must be adressed by two 
different lasers. One then chooses the laser resonant with 
(g; 2n, m) — > (e; 2n + 1, m) to propagate along +e y , and 
the laser resonant with (e; 2n + 1, m) — !• (g\ 2n + 2, m) to 
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propagate along — e y . According to Eq. (64), the sign of 



the phase factor changes on every second column, which 
results in a rectification of the staggered field to a uniform 
one (provided off-rcsonant terms are negligible compared 
to resonant ones). Practically, such a potential gradi- 
ent can be implemented by a real electric field affecting 
g and e equally, or by the AC Stark shift exerted by 
an off-resonant laser. The large electric fields or optical 
powers involved, combined with the fact that the tilting 
potential must be linear over the whole cloud to ensure a 
uniform transition freqency, make this option challenging 
from an experimental point of view. 

A possibly more practical configuration is based on a 
superlattice potential with period 2d x along x on the top 
of the main lattice, which can perform the same role as 
the potential gradient of the previous paragraph ( |Gerbier 



and Dalibard 20101. One needs three (instead of two) 



different coupling laser frequencies in this case, but this is 
not a significant increase in difficulty since the frequency 
differences are in the tens of kHz range, and can thus be 
easily gener ated using freq uency modulators. Another 
proposal by Mueller ( 2004 ) realizes the same goal by us- 



ing three internal states in three different sublattices. In 
that case, the distinction between neighboring transitions 
is automatic, and there is no need for an additional exter- 
nal potential. Unfortunately this trapping configuration 
seems difficult to implement in practice for the commonly 
used atomic species. 



E. Connection with dressed state approach 

It is interesting to establish a link between the opti- 
cal lattice case and the schemes discussed in sections U 
and[nj where we assumed that the atoms were following 
adiabatically one internal dressed state, defined as one of 
the local eigenstates of the atom-laser coupling matrix. 
We first discuss the case of lattices with laser assisted 
tunneling that we have considered in the first part of the 
present section. We then briefly p resent the con cept of 
flux lattices recently introduced by Cooper (2011). 



obtain the validity condition for the adiabatic approx- 
imation k (0) > y/V E R /h. In the TB limit the right 
hand-side is approximately the frequency gap w gap be- 
tween the ground (77 = 0) and the first excited (77 = 1) 
band. With the adiabaticity condition now written as 



,(0) 



> ui. 



g"I>- 



we immediately conclude that the adiabatic 
following of a given dressed state is intrinsically incom- 
patible with the single-band model used this section. In 
other words the laser-assisted tunneling scheme consid- 
ered in this section is a "diabatic process" Smith ( 1969 1] 
with respect to the dressed state basis, and it is therefore 
qualitatively different from the adiabatic scheme consid- 
ered in sections I and II to generate artificial gauge fields 
in a bulk system. 



Flux lattices [ Cooper] 2011). A flux lattice is based on 
a purely periodic 2D pattern of interfering laser fields. 
Atoms are modelled as two-level systems like in section 
|IJ with a coupling matrix U in Eq. ^ whose coefficients 
D, 6 and <j> are periodic functions of x and y. The 
atom-laser coupling is chosen large enough so that one 
approaches the regime where the atom follows adiabati- 
cally one dressed state (Eq. ( 17)), which is thus opposite 



to the case of laser assisted tunnelling discussed above. 
Because the vector potential A calculated from Eq. Q 
is also a periodic function of x and y, one could think 
that the flux of the magnetic field across a unit cell of 
the lattice, which is equal to the circulation of A on the 
contour of the cell, is always zero. However this does not 
hold if the laser field is chosen such that A has singu- 
larities inside the cell. Such singularities occur at loca- 
tions where sin# = 0, so that <j> is ill-defined. Contrary 
to A, the magnetic field B in these points, calculated 
from Eq. @, is non-singular and iCooperl rt201l[) showec 
that there exist configurations where B z keeps a constant 
sign over the cell, 2 with a magnitude ~ hk 2 . In particu- 
lar the Chern number of the lowest allowed energy band 



(Thouless et al. 1982 1 can be non-zero. This scheme thus 



constitutes an attractive alternative to the ones based on 
laser assisted tunnelling. 



Lattices with laser assisted tunnelling. The dressed states 
can of course also be defined in the lattice case. For 
simplicity we take a uniform Rabi frequency n^ ' for the 
coupling laser and a spatially varying detuning A(r) — 
Ao + (V e (r) — V g (r))/H, where V g and V e denote the po- 
tential energies for the g and e sublattices, respectively. 
We denote Vq the characteristic amplitude of oscillation 
of V e and V g , and we assume that fiK^ ' <C Vq to en- 
sure that the coupling laser does not lead to a strong 
deformation of the state-dependent lattice potential. As 
explained in Sec. |I.C[ this corresponds to a situation such 
that the mixing angle 9 varies rapidly around the points 
where the resonance condition A(r) = occurs. We use 
as a typical value |VA| m kV /H, where A: is a typical op- 
tical wave number. From the discussion of Sec. II. C[ we 



F. Non Abelian gauge fields in a lattice 

Optical lattices are also well suited for the generation 
of artificial non- Abelian magnetic fields, as first proposed 
by Osterloh et al. ( 2005 1 , and we conclude this section 



2 More precisely the total field consists of a regular background 
component with non-zero average, plus a periodic array of gauge- 
dependent Dirac strings — located in points where cos£> = — 1 for 
the gauge choice of Eq. (J3J - and carrying a flux of opposite 
sign as compared to the background. The regular and singular 
components together give a zero net flux over an elementary 
cell, as expected for a periodic Hamiltonian. However the Dirac 
strings are non-measurable and only the background component 
has an observable effect. 
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by a brief description of a possible implementation. As 
in Sec. |III[ the basic idea to generate non-Abelian poten- 
tials is to use several internal atomic levels that are cou- 
pled by various laser beams. Here we consider an atomic 
species with 2N quasi-degenerate sublevels, trapped ei- 
ther in sublattice g (sublevels gi, i = 1, • • • ,N) or in 
sublattice e (sublevels Sj, j = 1, ■ ■ • ,N). Instead of a 
simple phase, laser-induced couplings must induce a ro- 
tation in internal space, generated by a N x N matrix M 
with non-commuting components M x , M y , M z . Consider 
for instance a species with spin 1/2 in both the ground 
and excited manifolds (the spin 1/2 plays the role of a 
fictitious "color" charge here), which moves in an op- 
tical lattice that is now state-dependent both along x 
and y. Laser-induced tunneling is used along both axes, 
with the additional possibility of changing the spin index 
m z = ±1/2. For instance, the lasers inducing tunneling 
along y can be chosen to flip the spin (rn z — > m' z = —m z ), 
so that My oc aa x + /3a y , where a, f3 are numerical coef- 
ficients depending on the laser phases. The lasers induc- 
ing tunneling along x are then chosen to preserve the spin 
(m z — > m' z — m z ), so that M x cx jl + 5a z . Provided that 
<5^0, i.e. the x-lasers apply a different phase condition- 
ally on the internal "color", this realizes a non-Abelian 
potential. The physical effects that one expects from 
this non-Abelian settings are similar to those discussed 
in Section |III[ in particular the emergence of a Rashba- 



like spin-orbit coupling ( 


Dudarev et al. 


2004 


Goldman 


2007H2009||Goldman et ai.||2009a|b||Osterloh et aL[|2005 


Satija et al. 


poos]). 



V. OUTLOOK 

We have presented in this Colloquium the physical 
principles that lead to the generation of artificial gauge 
potentials on neutral atoms, using their coupling to laser 
fields. We have considered the cases of bulk systems and 
discrete lattices, and we have shown that both Abelian 
and non Abelian gauge potentials are accessible, provided 
one chooses a suitable atomic level structure and proper 
incident light fields. We have also explored several phys- 
ical consequences of these gauge fields, such as the nucle- 
ation of quantized vortices in a superfluid, the generation 
of a spin-orbit coupling via a non-Abelian gauge field, or 
the possibility to address the strong magnetic field limit 
in a 2D lattice with the corresponding Hofstadter but- 
terfly energy spectrum. The variety of the items in this 
(incomplete) list shows the richness of the situations that 
can be realized. 

These artificial fields also constitute novel tools for 
characterizing the properties of an assembly of atoms. 
Cooper and Hadzibabic| ( |2010[ ) suggested to use a small 
artificial magnetic field to probe the superfluidity of a 
gas. The field is generated with a laser scheme close to 
the one of Sec. |II.B[ and it is used to simulate a rotating 
bucket experiment and measure the reduced moment of 
inertia of a fluid when it acquires a superfluid compo- 



nent. In presence of the artificial field the normal com- 
ponent will stay at rest in the lab frame, whereas the 
superfluid component will rotate. A key feature of this 
proposal is that one can use spectroscopic methods to 
access the respective populations of the superfluid and 
normal states, by measuring the populations of the vari- 
ous ground states gj involved in the process. 

In order to keep this Colloquium within reasonable 
limits, we intentionally restricted our analysis to single 
particle or mean-field physics. However it is clear that 
these gauge fields can lead to very interesting phenom- 
ena when combined with the richness of strongly corre- 
lated states that emerge in many-body physics. For the 
Abelian case in the continuum limit, the situation is sim- 
ilar to the case of rotating gases, and paths to Quantum 
Hall physics that have been identified for Bose and Fermi 
gases remain of order [for a review see |Cooper|fl2008l )]. 
In the context of the Bose Hubbard model in presence of 
a uniform magnetic field, jMoller and Cooper ( 2009[ ) have 
recently proposed an approach based on the composite 
fermion theory to establish the existence of strongly cor- 
related phases that have no equivalent in the continuum 

Another line 



limit (see also |Polak and Kope6 



(20091) 



of research with optical lattices deals with the combina- 
tion of artificial gauge fields and nearest-neighbor inter- 
actions, as recently explored by Ruostekoski (2009 1 . 

To end on a somehow futuristic tone, we note that 
with the generation of non-Abelian gauge fields on an 
atomic gas, one will have at hand matter with intriguing 
topological characteristics. A possible application is the 
simulation of topological insulators with neutral atoms 
[see e.g 



Wu (2008), Stanescu et al. (2010) and Gold- 



man et al. ( 2010 )J . Another intrigui ng perspective could 



be topological quantum computing (|Sarma et al. 2006). 



In this respect the artificial gauge fields provide a differ- 
ent scenario compared to the anticipated emergent non- 
Abelian excitations in for instance fractional quantum 
Hall systems. With artificially created gauge potentials 
the pseudo-spin is effectively turned into a non-Abelian 
anyon, since two pseudo-spins which swap places will be 
represented by different final states depending on wether 
they were interchanged clockwise or counter clockwise. 
Such effects may provide an atomic building block for 
fault tolerant topological quantum computing. 
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